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Abstract 



^ , The Bakamjian-Thomas relativistic quark model provides a Poincare rep- 

resentation of bound states with a fixed number of constituents and, in the 

C^ , 

^SJ , heavy quark hmit, form factors of currents satisfy covariance and Isgur-Wise 

scahng. We compute the Light Cone Distribution Amplitudes of B mesons 
Qh' ^±{^) '^^11 as the Shape Function S{uj), that enters in the decay B — > ^^7, 

that are also covariant in this class of models. The LCDA and the SF are re- 
lated through the quark model wave function. The former satisfy, in the 



limit of vanishing constituent light quark mass, the integral relation given 
by QCD in the valence sector of Fock space. Using a gaussian wave func- 



' tion, the obtained S{u)) is identical to the so-called Roman Shape Function. 

I From the parameters for the latter that fit the B Xg'j spectrum we pre- 

. diet the behaviour of (/^^(a;). We discuss the important role played by the 

' constituent light quark mass. In particular, although (p^{0) 7^ for van- 

ishing light quark mass, a non- vanishing mass implies the unfamiliar result 
^ ■ </'-(0) = 0. Moreover, we incorporate the short distance behaviour of QCD 

' to ip^^uj), which has sizeable effects at large uj. We obtain the values for the 

parameters A = 0.35 GeV and A^^ = 1.43 GeV^^. We compare with other 
theoretical approaches and illustrate the great variety of models found in the 
literature for the functions 99^(0;); hence the necessity of imposing further 
constraints as in the present paper. We briefiy review also the different phe- 
nomena that are sensitive to the LCDA. The value that we find for X]^^ fulfills 
the upper bound recently measured by BaBar. 
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1 Introduction. 



The Light Cone Distribution Amphtudes (LCDA) of heavy-hght mesons ^^{C,), 
$:^(^) [1] or, in the heavy quark hmit if^{uj) and f^{uj) [2j, are fundamental func- 
tions that enter in the large energy limit of amplitudes of semileptonic decays [2] 
and in non-leptonic decays of B mesons [3], in the determination of the form factor 

^^(0) [1] and, more directly, in the decay -fiui |5l El [3 E] . 

On the theoretical side, these functions have been studied in the perturbative 
regime at large [Til^, and a number of very varied Ansatze have been proposed for 
the dominant part of them at low ^ or u, where the function is peaked at ^ ~ ^^^^ 

orcu- Aqcd [21 mi- 

On the other hand, one can obtain model independent information on these 
functions from the measurement of the spectrum in the decay B^ 'yi'J^e, that is 
directly related to one of the LCDA [3 El [3 [8] . 

On the theoretical side, although rigorous results are known in the perturbative 
regime uj ^ ^qcdi the guesses advanced for the main non-perturbative part of 
the LCDA amplitudes come essentially from QCD Sum Rules, imposing continuity 
between the perturbative and long distance regimes [TJ [9] . 

The motivation of the present work is to compute the LCDA in a class of rela- 
tivistic quark models, namely the Bakamjian-Thomas (BT) quark models [101 ttH 
W]\ |T31 HI] . This is a class of models with a fixed number of constituents where the 
states are covariant under the Poincare group. The model relies on an appropriate 
Lorentz boost of eigenfunctions of a Hamiltonian describing the spectrum at rest. 
On the other hand, one has demonstrated that the matrix elements of currents be- 
tween hadrons are covariant in the heavy quark limit and exhibit Isgur-Wise scaling 
[T5] in this limit [H] . Given a Hamiltonian describing the spectrum, the model pro- 
vides an unambiguous result for the elastic Isgur-Wise function i{w) jTH [16]. On 
the other hand, the sum rules (SR) in the heavy quark limit of QCD, like Bjorken 
or Uraltsev SR are satisfied in the model [TTl |18], as well as SR involving higher 
derivatives of ^(u;) at zero recoil [T9] . 

The interest of computing the LCDA functions in this framework is to 

directly relate them, in an unambiguous way, to the Shape Function Siu) [201 EH 
[221 [23l [2^1 [251 [26] 'm. B ^ ^s7, that can also be computed in the BT class of models. 
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One could calculate the LCD A within the BT scheme using a Hamiltonian de- 
scribing the spectrum, like the Godfrey-Isgur (GI) model [27j, that has been used 
elsewhere to compute the elastic Isgur-Wise function C,{w) and the inelastic ones 
Ti/2{w), T3/2{u)) [16]. However, we have tested the GI model to compute the Shape 
Function S{uj) and have realized that this model does not fit the B Xg'y spec- 
trum. This is the reason why we have decided another phenomenological approach, 
namely to relate the Light Cone Distribution Amplitudes ip^{uj) to the Shape Func- 
tion S{uj), a relation that is provided by the BT scheme. Using this relation, and a 
SF ^(ti;) that fits the B — > Xg'j spectrum, one can predict the LCDA. The problems 
for the SF and also a discussion of heavy-to-light form factors within the BT scheme 
for the GI spectroscopic model will be given in detail elsewhere. 

The paper is organized as follows. In Section 2 we give the master formulae 
defining the theoretical framework of BT quark models. In Section 3 we review 
the definitions of the LCDA at finite mass and in the heavy quark limit of QCD 
and in Section 4 we set a main natural hypothesis to compute the LCDA within 
quark models. In Section 5, to introduce the technicalities of the BT model, we 
review the calculation of the heavy meson decay constant. In Section 6 we obtain 
our main results, namely the expressions for the LCDA in the BT quark models. In 
Section 7 we compute the SF ^(ci;) in BT models and show that, in the case of the 
harmonic oscillator, it is identical with the so-called Roman Shape function, used 
to fit the B — > Xs7 spectrum [Sni EH [22l [26] . In Section 8 we use the parameters of 
the latter to predict the LCDA functions (p±{uj) and their moments. In Section 9, 
following Braun, Ivanov and Korchemsky [7J and Lee and Neubert [9] , we introduce 
the radiative tail of (p^{uj). In Section 10 we compare our results with proposals 
for the LCDA in other theoretical schemes. In Section 11 we review the different 
phenomena that are sensitive to the LCDA, and in Section 12 we conclude. 



2 The Bakamjian-Thomas relativistic quark model. 

As explained in [H] , the construction of the BT wave function in motion involves 
a wnztar?/ transformation that relates the wave function s„(Pi5 ' ' ' ' P™) terms 
of one-particle variables, the spin Sj and momenta Pi to the so-called internal wave 
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function ^^(P, k2, • • • , k„) given in terms of another set of variables, the total 
momentum P and the internal momenta ki,k2, • • • ,k„ (X^kj = 0). This property 

i 

ensures that, starting from an orthonormal set of internal wave functions, one gets an 
orthonormal set of wave functions in any frame. The base 1.. ^^(pi, ■ • • , Pn) is use- 
ful to compute one-particle matrix elements like current one-quark matrix elements, 
while the second \E'*'^*... ^^(P, k2, ■ ■ • , k„) allows to exhibit Poincare covariance. In 
order to satisfy the Poincare commutators, the unique requirement is that the mass 
operator M, i.e. the Hamiltonian describing the spectrum at rest, should depend 
only on the internal variables and be rotational invariant, i.e. M must commute 
with P, ^ and S. The internal wave function at rest (27r)^(5(P)</?s^^...^s„(k2, • • • , k„) 
is an eigenstate of M, P (with P = 0), and S;^, while the wave function in mo- 
tion of momentum P is obtained by applying the boost Bp, where P° = VP^ + M"^ 
involves the dynamical operator M. 

The final output of the formalism that gives the total wave function in motion 
^sr,- -,s„(Pi' ■ ■ ■ ' P") terms of the internal wave function at rest (ps^,...,sn{}^2i ■ ■ ■ , k„) 
is the formula 



*s..,„(Pi, • • • ,Pn) = {2nrs f ep. - p 1 ^ I n 



where p° = y^pf -|- m| and Mq is the free mass operator is given by Mq = ^{Y.PiY- 
The internal momenta of the hadron at rest are given in terms of the momenta 
of the hadron in motion by the free boost ki = Bv^ pi where the operator Bp is 

i 

the boost (VP^, 0) p, the Wigner rotations Rj in the preceding expression Rj = 
B-/B^^.Bfc, and the states are normalized by < P'^S'^lP^S^ > = {27rf5{P' - 

The current one-quark matrix element acting on quark 1 between two hadrons 
is then given by the expression 

*<'-,^n(P'l' • • • 'P")* < p'l,s'i\J^'^\pi,Si > (Pl, • • • ,Pn) (2) 
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where ... s„(pi, ■ ■ ■ , Pn) is given in terms of the internal wave function by ([T]) and 
< p\, s'i| J^^^Ipi, Si > is the one-quark current matrix element. 

As demonstrated in [Ml [28], in this formalism, in the heavy quark limit, current 
matrix elements are covariant and exhibit Isgur-Wise scaling, and one can compute 
Isgur-Wise functions like ^{w), ti/2{w), T3/2{w) [16]. 

In the present paper, as far as the LCDA are concerned, we are dealing with 
current matrix elements between one meson and the vacuum, i.e. < 0|J|\I' >. In 
[29] such matrix elements were considered and it was demonstrated that the decay 
constants of heavy-light mesons are covariant - independent of the frame - in the 
heavy quark limit, and exhibit heavy quark scaling, i.e. /sy^me is a constant in this 
limit. This quantity was calculated for various Hamiltonians describing the meson 
spectrum. 

We want now to go beyond and compute the LCDA '^±{uj) starting from the 
meson-to- vacuum matrix elements < 0| J|\I^ >. We can here advance our main result, 
that is parallel to the one obtained for meson-to-meson current matrix elements. The 
LCDA are covariant in the heavy quark limit, and can therefore be computed without 
any arbitrary parameter once the Hamiltonian giving the internal wave function is 
given. The same statement holds for the Shape function S{uj) in B Xgj, that can 
be expressed as a meson-to- meson matrix element < \l/ 1 1 ^ > . For the latter we will 
use the harmonic oscillator that gives the Roman Shape Function [20l [2T1 [22l [26] . 
We will then use the phenomenological parameters of this function, that fit the 
B — > spectrum, to predict the LCDA. 



3 B meson Light Cone Distribution Amplitudes 

Let us define the LCDA $^(0, [D 

<0|g(2)5„_(z,0)r6(0)|5,(P) > 



V ■ n_ 



(3) 



where 5'„_ (z, 0) is the Wilson line following the light-like four vector n_ = (1, 0, 0, —1) 
(n^ = 0), f is the B four velocity v = F is an arbitrary Dirac matrix, and the 



center-of-mass motion is along the Oz axis. P+, z_ are light-cone variables defined 
for any four vector p by p+ = ^ , V- = ■ Sometimes one takes v ■ = 1, 
i.e. the -B rest frame, but to exhibit covariance we have adopted a general value for 
V -ri-. The LCDA (^bAC), '^'^2(0 satisfy the normalization conditions 

('di^l{i)= fdi^^_{0 = l (4) 
Jo Jo 

In the heavy quark limit — > oo it is useful to use a new variable 

uj = rub^ (5) 
keeping u fixed, that yields the definition of the LCDA [21 HJ [30] . 



fsmB 



< 0\q{z)S^_{z,0)TK{0)\Bd{v) > _ (6) 



1^ du e-'"''+"-Tr |r(l +y)75 



4 Jo 

The relation between $^(0 v'^(co') is 



v^+M 



V ■ n_ 



^1(0.) = (^) (7) 

and v^^(c(;) satisfy the normalization conditions 

du (^l(cj) = 1 (8) 



oo 







4 LCDA in quark models. 

In what follows, we will use the preceding relations to obtain the expression of 
(/?^(cc;) in BT quark models. In the class of BT quark models, gluon exchange is 
included in the potential. Therefore, in a way, the gluon field is integrated out, but 
one looses the explicit gauge invariance that is ensured by the Wilson line of the 
preceding expressions. 

In a quark model, what we can consider is the matrix element involving con- 
stituent quarks, in particular constituent light quarks with a dynamical mass. Our 
ansatz will be to identify the QCD matrix element with the Wilson line with a ma- 
trix element involving the constituent quark field, or in more rigorous terms, one 
would say that one works in the light-cone gauge, 

^+ = 5„_(;2,0) = 1 (9) 



6 



and set 



< 0\q{z)Sn_{^,0)rb{0)\Bd{v) > ^^^^^^^ 
= < 0|g,_,,„(z)r6(0)|5,(t;) > (10) 

This is our main hypothesis and the starting point of the quark model calculation, 
that then follows in a straightforward way. From now on the constituent light quark 
field qconstituent will be denoted by q. Of course, the condition ^ can hold only in 
field theory, and our BT scheme is just a model. Defining $±(0 by the quark model 
expression : 

<0|g(^)r6(0)|S,(P)> 



- Y e-^«^+-Tr |r(f + 771^)75 



V ■ n. 



where n_ = (1, 0, 0, —1), one obtains 



fsmB 27r 



J V ■ 



z+=zx_=0 



fBrriB 271 J \ vn_) 



(12) 



Calling P2 the four- momentum of the light quark, using translational invariance 
and integrating over z_, these expressions write 



1 



fBrriB 



< 0|g(0)<5 



P2+\ i- 



P< V ■ 



75K0)|5,(P) > 



<0|g(0)5k-^ ) 2y- 



V ■ n_ 



75&(0)|fid(P)> (13) 



These will be the starting formulas to compute these functions in the BT quark 
model, from which we will deduce their heavy quark limit !f^{uj). But let us first 
compute the heavy meson decay constant Jb in the BT quark models, that will 
provide the desired normalization for the LCDA. 



5 B decay constant in BT models. 

The calculation of the matrix elements to obtain the LCDA ^p^{uj) is just rem- 
iniscent of the one made to obtain the corresponding decay constant [29j. In this 



latter case one needs the matrix element 



<0|g(0)r6(0)|5,(t;)> 



dp2 EiPi 



Tr 



k^k2 jmim2 1 



<^(k2) 



fci u 



2 " 2 

where </5(k2) is the internal wave function at rest, with the normalization 

dk2 



1 



(14) 



(15) 



Pi and p2 {mi and 1712) are the quark four-momenta (masses ) of respectively the 
heavy and light quarks, Bp is the 4x4 boost matrix associated with the four- vector 
p, and the four vector u, Mq and the relation between ki and pi are given by 

Pi +P2 



u 



Mn 



Mn 



(Pl +P2Y 



^uki =Pi{i = 1,2) 



(16) 



where ki and ^2 are the four-momenta of the quarks in the rest frame of the B 
meson and B„ is the boost associated with the four-vector u. The products of 4 x 4 
matrices under the trace read 



B„Bfc2 — ^ — h}„ 



m2 +^2 



2m2 (^2 + ^2) ^ 
l+'^i mi 



2mi {ki + mi) 



(17) 



This yields the expression 



<0|g(0)r6(0)|5rf(P)> 



(2vr) 



1,0 /,0 
1^1 1^2 



p\pI\ {kl + mi) (A;^ + m2) 



Tr 



r'T (m2 +^/2) (1 +i) (mi (/?(k2) 



Using the current V = 7^75 we obtain, after some algebra and the definition of the 
four vector u ffT6l). 



<0|g(0)7^75&(0)|i?d(P) >= 



1,01,0 
"'1"'2 



v/2 i (27r)3pO|,o\ (fco + mi) (A;2° + m2) 



mi + m2\ 
1 + — ) (Pi/.m2 + P2i,mi) ip{k2) 



(19) 



Since the BT states are normalized acccording to 



< Ba{P')\Bd{P) >BT = (27r)3(5(P - P') (20) 

while the covariant normalization is < 'Bd{P')\Bd{P) >= (27r)^2P°5(P - P'), we 
have to identify the former matrix element with the definition of the decay constant 

< 0\q{Oh,^m\MP) >= ^ (21) 

one obtains 

dp2 Vv^ 



/C°/C2 f^^mi+ 1712 



(27r)3p0p0^ (A;0 + mi)(A;0 + m2) V Mq 

[mi (P2 -V) +1712 {pi ■ v)] Lp(k2) (22) 

This expression is not covariant, but becomes covariant in the heavy quark limit. 
For mi ^ oo one has 

Pi 

mo mi ^ niB , u v , > v , — > 1 , kn ^ P2 ■ v (23) 

mi mi 

where v is the B meson four-velocity, and one gets the expression of the B decay 
constant in the BT model in the heavy quark limit |29j . 

JbV^ = V2v^ / ■^^^^JiP2-v){p2-v + m2) ^ [^iP2-vf-ml^ (24) 

This expression is covariant, satisfies heavy quark scaling and gives the decay con- 
stant in terms of the internal wave function. In the B rest frame one gets 

1/2 



We have checked that these expressions for the decay constant hold exactly in 
the equivalent light-front approach of Cardarelli at al. [13] in the heavy mass limit. 



6 LCDA in the heavy quark hmit in BT models. 



According to f[T3|) . we have to compute a generic matrix element 

= ^ < mm - ^1 m)mp) > (26) 
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2+ = Zj^=0 



where F is a Dirac matrix. In the BT quark model, this expression writes, taking 
into account the appropriate normahzations. 



^^2- 



fBTTlB 8 J (27r)3pO 

-Tr 



(27) 



1,01,0 



7T(m2+^/2)(l+^)(mi+yi) yp(k2 



\\ (A;? + mi) {k^ + ms) 
Changing the measure and integrating relatively to p2+ and p2-, one obtains 

1 ^/N'rl 1 



$(0 = -V2P0 



fBTUB 8 e27rLi (27r)2' \ (A;? + (A;0 + ma 



lOlO 



Tr 



7^r(m2+^2)(l+^)K+^^i) '^(^2 



(28) 



2P+i 



Making use of the definition of the four-vector u (fT6|) and computing the trace 
particularizing respectively to F = 75 and T = (2^— -^^^ 75, one gets 



mi + m2 
Mn 



1 \/iv:i 1 



/^ms 2 ^27rU (27r)2^ ^ (A;0 + (A;0 + ma) 



1,0 7,0 



mi{p2 ■ ^-) + "^2(^1 ■ n-) 



V ■ n. 



V5(k2 



(P2_L)^+m| 
P2+=5P+,P2- = 2P^.g 



^rl 1 



2± 



1,01,0 
/ti/t2 



/^m^ 2 ^ 27r 

^ mi + m2 



^ (fc? + mi) (A;^ + m2) V" ' Mq 

"^1(^2 ■^-) +"^2(Pi ■'^-)1 
) ka 



{2 [mi(p2 ■ f ) + m2(pi ■ v)] 



(P2X)^+" 

2P, e 



(29) 



6.1 Heavy quark limit. 

At finite mass the expressions ( l29i) are not covariant. In the heavy quark limit, 
using now the variable uj = mi^. 



U ^ Vi ^ V 

Mq nil i^B 
one obtains, denoting p2_L = p_l 



/c° — > mi , k2 ^ P2 ■ V 



(30) 
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2uj 



\ p]_ + {00 + 7712)2 



\J (pi + + rnl)^ — Aijj'^m'2 



(31) 



In the heavy quark hmit, the LCDA ^^{00) are covariant since the boost is along 
Oz, and the variable uj can be written in the covariant form, 



^ P2+ P2 -n^ p2- V. 

u = vfiBt, = rriB— = rriB— = 

r ■ U- V ■ ri- 



(32) 



Performing the change of variables 



P 



'pj_ + + ^2) — Auj'^ml 
2uj 



(33) 



one obtains 

B 



ip+{uj) = - — -p^V2—r{uj + m2) dp 



P 



Tj2'P(P) 



dp- 



p 2Jp^ + 1712 — UJ + m2 



with 



Po[^) 



(^2 + miy/'^ ^ _|_ 77^2 + 7r72 

Tj2Vip) (34) 

(35) 



(p2 + 7712)^/"^ [yP'^ +1712 + m2 
\iU - 77721(^0" + 7772) 



2UJ 



One can check the normalization ([8]) by changing the order of the integrals 



UJ > 0, p > Po{uj) 4^ P > 0, uj_{p) < UJ < uj+{p) 



with 



UJ±{p) = Jp"^ + 7172 i P 



(36) 



(37) 



The integrals over uj are trivial 

fUJ+ip) f^^+(p) 



I duj{m2 + uj)= I duj i^Jp"^ + 7172 — UJ + ^2) = 2p (^p'^ + ml + m2 



(38) 



and the normalization (IH]) follows. 
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6.2 Limit of vanishing light quark mass. 

In the limit of vanishing hght quark mass, one can immediately demonstrate from 
expressions (1341) that ^p±{u}) satisfy the differential equation, the so-called Wandzura- 
Wilczek approximation [1] [3T] 



or, equivalently, the integral relation 



(39) 



^«(.)=r,p^5M (40) 

Jtjj p 

Moreover, for vanishing light quark mass, one gets from for small uo, 

if^ico) ~ u ^^(uj) ~ Constant (m2 = 0) (41) 



The relation (1391) or ( 1401) holds in QCD if one restricts to the valence quark sector 
of Fock space |31 [31j. It is reassuring that this relation holds also in the BT quark 
model, since it is obviously formulated in the valence quark approximation. For this 
relation to be satisfied one needs nevertheless the dynamical light quark mass to 
vanish. However, it will be interesting to grasp the significance of the corrections 
due to a non- vanishing m2 and check its numerical effect on ip^{uj), as we will do 
below. 

6.3 Behaviour of ^±{uj) for non- vanishing Hght quark mass. 

A specific trend of our results for the LCDA ^p^{uj) is the important role played 
by the light quark mass m2, as will show the numerical results of Section 8. The light 
quark mass has dramatic consequences, namely that the first derivative of ^p^{uj) 
and f^{uj) vanish, as we can see by inspection of formulas ( 134|) . (I35l) : 

^f(0) = ^^(0) = (42) 

This is apparently at odds with the general belief, based on QCDSR, that predicts 
the behaviour ( l4T1) [21 [6l El [8] . The relation ( l40l) is strongly violated in the presence 
of a constituent light quark mass, since the l.h.s. of (HOl) vanishes : ^p^{uj) vanishes 
for UJ ^ 0, since the lower limit fl35|) of the integral behaves like 



2 

Po{uj) ~ > oo for a; — > (43) 

UJ 

12 



In the QCDSR approach one has an expression giving the correction to the 
behaviour (jUj) from the < qq > condensate, i.e. the first order correction due to 
the dynamical hght quark mass, the gap induced by dynamical chiral symmetry 
breaking. This is the constituent mass involved in our quark model calculation, 
proportional to the condensate. The corrections at first order in < gg > have been 
given by Grozin and Neubert [2] that write, in a simplified notation, 

^^{u) = ^^J'\co)- <qq>M (44) 

where ^±{uj) are LCDA at leading order, independent of < gg >, satisfying (1391) . 
and the function f{uj) has the behaviour 

7(0) = /(oo) = (45) 

Notice moreover that the correction dependent of < gg > is the same for ^^iyj) and 

In our scheme, the first order corrections in 1112 to and {uo) are also equal, 

but our calculations contain higher orders in the constituent light quark mass, as 
we observe from fl3Tl) . flMl) . fl35|) . giving the much stronger behaviour P2l) . 



6.4 Moments of ip^{uj). 
Defining the moments 

Mf ) = / du u^'^liu) m1^) = / du cu^y.!(a;) (46) 

Jo Jo 

one finds for m2 7^ that M^^ = mI"-* = 1, i.e. the normalization condition ([8]). 
Moreover, one finds for any moment with > 0, for vanishing light quark mass 
7712 = 0, 

Mi^^ = {N+ l)Mi^) (47) 

that holds obviously in the BT class of quark models since it follows from the QCD 
relation fl39|) . Of interest are the moments M^^ that are given in the valence 

sector of QCD by [2] _ 

M« = 2M^y = — (48) 
3 
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and should allow to compute A in the BT models. 

The moment = A^^, satisfies, due to the positivity condition in the absence 

of a radiative tail (p+{uj) > |5], 

7 The Shape Function S(ct;) in the BT approach. 

7.1 Definitions of the Shape Function. 

The shape function S{uj) that enters in the description of the decay B — > Xgj 
[22l [23] is defined by the expression 

S( ) 1 Tr (t^) = < B{PB)\hvr6{uj - zD+)K\B{pb) > 
^'^^ 2 V 2 y 2mB 

where p+ = p° and the support of S{uj) is 

-oo<cj<A (51) 
One uses also another definition [231 [25] 

^(cD) = S{A - Co) (52) 

with the support 

< cD < cx) (53) 
The functions S{uj) or 5'(cD) can also be computed in the BT model. 

7.2 Calculation of the Shape Function in the BT model. 

Since the field annihilates the b quark within the B meson, one can introduce 
a complete set of intermediate states of the spectator quark \p2, S2 >< P2, in the 
preceding expression 

< B{pb)\K{x)6{uj - iD+)K{x)\B{pB) > = 

<B{pB)\K{x)\p2,S2><P2,S2\S{uj-iD+)K{x)\B{pB)> (54) 

P2,S2 
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with 

<P2,S2\K{x)\B{pb) >=e'^P'-P^+'^'^>^ <P2,S2\K{0)\B{pb) > 

< B{pB)\Kix)\p2,S2 >= e-'^^^-P^+^^^y-' < BipB)\KmP2,S2 > (55) 

We use now for the operator 6{uj — i-D+) the following representation 

1 

6{uj - iD+) = — ds e'("-'^+)^ (56) 

2n J-oo 

and make the identification in the quark model, that follows from the hypothesis (Q 

D.^O^-i, (57) 
One obtains, for any function f{x~^) 

OO 1 

from p- X = p+X- — p± ■ x_l with p± = p^ ±p^ and 9+ = = one finds 

gs9+ ^^{p2-PB+m^,v)■x 

— g*{(P2-Ps+mi,t)) + (a;_+s)+{p2-Ps+m6i))-a;+-(p2-ps+mbv)x-Xx} |-gg^ 

and therefore 

5{UJ - ^5+)e^(P^-^^+'"''">" =5[UJ+{P2-PB+ nH,v)+] e^V2-VB^m,v)-x ^gQ) 

hence, in the quark model 

S{u) - Tr fr-^ + - ^ '^^VB)'^v{x)b{uj-%d^)K[x)\B{j)B) > 



2 \ 2 J 2mB 

< B{pb)\K{x)5 [uj + {p2 -PB + mhv)+] K{x)\B{pb) > 



(61) 



2mB 

In the B rest fame mbV+ = rrih, ps = {itlb,0) and therefore, from A = rriB — mb one 
gets 

S(u;) - Tr (r^) = < BipB)\Kix)r5{uj + P2+ - A)K{x)\B{pb) > 
^ ^ 2 I 2 1 2mB 



One obtains for the shape function 

S(Co) - Tr (r^) = < BipB)\M^)T6{Q-p2+)Kix)\BipB) > ^^^^ 
2 \ 2 I 2mB 
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In these expressions p2+ is the component p2+ = P2+P2 of ^^e spectator hght quark. 

Let us now compute exphcitly this last expression in the BT model, proceeding 
along the same lines as we have done in Section 6 for the LCDA. One needs to 
compute the matrix element of the operator r6{Q — p2+) in the forward direction 

dp2 1 



<fi(P)|r5(cD-p2+)|5(P) >= 



(27r)3 



[(^(kz)] 6{uj-p2+] 



(64) 



u 



01,0 
2 



-Tr [r(mi +i){m2 +^)(mi 



(p?) (fc? + mi)(fc0 + m2) 16' 
where ki are related to pi through the boost ( |T6l) . B^/cj = pj. Changing the measure, 
and taking into account that now p2+ = P2+ P2-> can write 



< B{V)\T5{uj-p2+\B{V) >= 

dp2+ dp2- . / p5_ + TTin \ dp± M° 

5 P2- - 



P2+ 27r 



P2H 



(27r)2 (pO) (A;? + mi) (A;^ + ms) 



(65) 



[y.(k2)]' 5(cD - P2+)^Tr [Tim, +^^(1 +^)(m2 +^)(mi 
Integrating relatively to p2+ and p2- 

< 5(P|r5(cD-j92+)|5(P) >= 1 -1 ^ ^ 



16 27r7 (27r)2 Q 



b(k2 



lOlO 
"'1"'2 



^ (rf)2 (A;0 + mi)(A;0 + m2) 
Tr [r(mi +yj(l +^)(m2 +^/2)(l +^)(mi 



(66) 



P2+=UJ,P2- 



7.3 Heavy quark limit. 

In the heavy mass limit one has u ^ vi ^ v , k2 —>■ P2 ■ v and Mq 
and therefore, after some algebra, 

<5(P|r5(cD-p2+)|5(P)>= 1 1 Mp± 1 



8 27ry (27r)2 tD 
P2 ■ 



P2 ■ V + m2 



Tr [r(l+y)(m2+^/2)(l+y)]} 



P2+=U),P2- 



and from 



Tr [r(l +y)(m2 +^2)(1 = 8 K + {P2 ■ v)] ^Tr 



mi — > mB 



(67) 



(68) 



16 



one gets finally 



<»-..,P(P>4r.(ri±^)±/|^i 



LP U {j)2 ■ vf - ml 



-XP2-v) 



P2+=UJ,P2- 



(69) 



Let us identify with the definition (152!) . With the normalisation of the BT model 
( !20|) we have to identify the matrix element (!69|) with 



2Mb 2 ^\ 2 y 2P0 



(70) 



Therefore one gets 

^ 1 f dp± 1 

and from 



2 7 (27r)2 uj 



{P2 ■ V) 



P2 ■ V 



p5_ + ml + tD^ 
2Q 



(71) 



(72) 



one obtains finally 



1 r dp_L p5_ + m| + tD^ 

2 7 (27r)2 25 



25 



(73) 



Again, similarly to the LCDA, one can obtain a more compact form of S{Cd) by 
peforming the change of variables 



k 



20 



(74) 



and one obtains the simple final result for the Shape Function in the BT model 



SiQ) 



with 



ko(Q) 



dk k[ip{k)Y {0<Q <oo) 



\Q - m2\{Q + m2) 
2Q 



that gives the value cDmax for which S{u)) attains its maximum value S„ 



^max — ^2 



47r2 Jo 



dk k[ip{k)Y 



(75) 
(76) 

(77) 
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For the other definition of the shape function S^u) (l50l) one reads, from (152!) . 



1 

47r^ JKo{io) 



with 



= I -^i^- ^ (-00 < ^ < A) (79) 

To check the normahzation condition 

du) S{Q) = 1 (80) 

we use ( !78|) and exchange the order of the integrals over cD and k, 

tD > , k > ko{uj) <^ A; > , Q^{k) < Co < Q+{k) (81) 

where 



Q^{k) = Jk^ + ml±k (82) 



One obtains 



^ 1 r^+(^) 1 r°° 

. dQS{Q) = — dkk[ip{k)f dQ = — dk P[ifi{k)]^ = 1 (83) 

47r^ JO J^^(k) 2t{^ Jo 

from the normahzation of the internal wave function (fT5|) . 
From (183!) and (!52|) it follows, for the function S'(co'), 

^ du S{uj) = 1 (84) 



7.4 Gaussian wave function. 

In principle, one could calculate the wave function f{k) entering in ip^{uj) (!34|) 
and S{uj) (!78|) from the quark potential. However, as argued in the Introduction, we 
prefer here to adopt simply a gaussian wave function, that has been used elsewhere 
to obtain S{uj) and to fit the B X,7 spectrum [2011211 1221126] . 

Let us compute the shape function S{uj) using the harmonic oscillator potential, 
with the internal wave function at rest v^(k2) is given by 

^Ho(k,) = (2vr)3/2 (v)'^'^^P (-^) ^85) 
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We obtain 



R 



exp 



(A - - ml 



4 i?2(A-cu) 



With the notation 



Pf 



Pf 



one gets 



1 1 



exp 



1 
4 



Pf P A ^ 

^ -^ (1 ~ ^) 

A 1-x 



(86) 



^7) 



that is correctly normahzed to 1. 

We obtain therefore in the BT model with harmonic oscillator potential the 
so-called "Roman" Shape Function [201 EH ESI EE] . 



8 Predictions for the LCD A (p±{uj). 



8.1 Parameters of the Roman Shape Function. 



Limosani and Nozaki [26] have made a recent fit to the Belle B Xs'j data 
[32], in order to extract the 6-quark shape function parameters Asf and Xf^ {Xf^ = 
— /i^) in a number of models for the shape function, among them the Roman Shape 
Function. Writing fl87p in their notation 

00 = k+ F{k+) = S{uj) (89) 



F(M = iV^exp|-l 



K 1 — X 



nil — x) 



(90) 



with 



0; = ^ . = ^ p = ^ (91) 

Asf Pf Pf 

amounts to replace in flSSl) A hy Asf of the Shape Function renormalization scheme 
[33]. Of course, the BT quark model is not field theory, and this scheme is too 
rough to distinguish between A and Asf, that differ by QCD radiative corrections. 
Therefore, we can make the replacement of A by Agi;' in the formulas of the precedent 
section, and this will be our final model for the shape function. 
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The first moments of F{k^) 

An= dk^ kl F{k+) 

J —oo 

are given in QCD by 

\SF 

Ao = l = A2 = 1- 

Computing tliem explicitly using the shape function ( l90l) one obtains 



k 

and the conditions (^5]) yield respectively 

1 



N 



A 



SF 



2 (1+P 



(Asf)' 



The best fit of Limosani and Nozaki |26] gives 

A^^ = 0.66 GeV Af^ = -0.39 GeV^ 

that corresponds, using (jHSD, to the values 

p = 0.776 K = 1.462 

or, from (IHTI) and (!9T]) . to the quark model parameters 

= — = 2.216 GeV"^ ms = 0.398 GeV 

Pf 

On can also compute the value (fc+)maa; of k+, for which F{k^) becomes 

{k+)max = ( 1 ^ \gp Fmax = N —= 

\ K J y/TT 

or numerically, 

ik+)max = 0.262 F^ax = 1-250 GeV-1 
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The other definition of the Roman shape function S{Q) (l52l) has a very simple 
expression in the BT model 

\ 2' 



S{Q) 



R 



exp 



LU 



that is identical to S{(2i] 



Roman 



with 



and 



p = R^ml 



A 



SF 



exp 



SF 



0.398 GeV 



Sm^a. = 1.250 GeV-^ 



(101) 



(102) 



(103) 



8.2 Predictions for (p^{u). 

Using the harmonic oscillator wave function (1851) . the parameters ( I98l) and ex- 
pressions (|3T|) or flMl) we can predict the LCDA ^p^{uj) within the BT model. 

In the limit of vanishing light quark mass 1712 = 0, in which the relation fl39|) 
between f^{uj) and v^:?(co') holds, one obtains simple analytic expressions in terms 
of the Error Function $(x) 

2 



vr Jo 



dt e 



(104) 



For the decay constant one gets 



fB^/rn^ 
and for the LCDA, 



— / du u 



TT 



1 IR 



3/4 V 2 



1 - $ 



(105) 



/bv^ vr3/4 V 2 



1 - $ 



(106) 



'2 / R^uj^' 
exp 



1 - $ 



To illustrate numerical results, and give a feeling of the dependence on the light 
quark mass, we chose two sets of parameters : 
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12 3 4 

Figure 1: The function ^p^{uj) for a gaussian wave function for the sets of 
parameters ( l98l) obtained from the Roman Shape Function [26] {R, m2) = 
(2.216 GeV-\ 0.398 GeV) (higher curve) and {R,m2) = (2.216 GeV~\ 0) 
(lower curve). For non- vanishing m2 the first derivative vanishes at a; = 0. 



(1) The reahstic case of the parameters of the quark wave function (|98l) R = 
2.216 GeV^^ and m2 = 0.398 GeV. From these parameters fitted on the B — > Xg'j 
spectrum we predict 

fBy/^ = 0.388 GeV3/2 (107) 

and the functions (p^{uj) and ip?-{uj) that are plotted respectively in Figs. 1 and 2. 
The heavy quark limit value (11071) gives, using the physical B mass, Jb = 170 MeV, 
a little smaller than the popular values for this quantity. 

(2) For comparison we adopt the vanishing light quark mass case, taking the 
same radius R = 2.216 GeV~^ but ms = 0. We find fs^/^ = 0.315 GeV^/^ ^nd 
the functions ^p^{uj) and ip^{uj) that are also plotted in Figs. 1 and 2. 

We observe that the non- vanishing light quark mass gives (p ^(0) = Oandy^f (0) = 
0, while for vanishing light quark mass one obtains (^^{0) ^ and also f^^' (0) 7^ 
like in QCD Sum rules. Since the vanishing (f^{0) = for m2 7^ is an unfamiliar 
feature, we plot in Fig. 3 the evolution of (p^{uj) with increasing values of 1712 = 0, 
0.1, 0.2, 0.3 GeV. 
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Figure 2: The function (p^{uj) for a gaussian wave function for the sets of parameters 
obtained from the Roman Shape Function [26] {R,m2) = (2.216 GeV^^, 0.398 GeV) 
and {R,m2) = (2.216 GeV~^, 0). For non-vanishing m2 the function vanishes at 
u = 0. 




Figure 3: Evolution of <f^{uj) with the constituent hght quark mass for R 
2.216 GeV-^ and ms = 0, 0.1, 0.2, 0.3 GeV. 
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8.3 Moments of (p^{u). 

Let us now compute the moments of ^p^{uj) using the harmonic oscillator wave 
function fl85|) that gives the Roman Shape Function fl90|) . fl9Tl) with the parameters 
( l98l) . We obtain for the first moments M^^ 

Mj^)^^ = 2M^y^'^ ^ = 0.964 GeV (108) 

3 

that gives Asf = 0.723 GeV, in qualitative agreement at the 10% level with the 
value fl96|) obtained from the fit to i? ^ Xg'y, and for the moment M]_ , called 
also A^^ 

[Xb^Y^ = M^'^^^^ = 1.521 GeV-^ (109) 



9 Radiative corrections. 

The importance of the radiative tail has been underlined by Braun, Ivanov and 
Korchemsky [7j. When taking into account one- loop QCD corrections, the moments 
for > of ^p^{uj) are divergent, as already pointed out by Grozin and Neubert ^ 
and by Descotes-Genon and Sachrajda [6J. 

Braun et al. give a parametrization of ip^{uj,fi) based on QCD Sum Rules plus 
the QCD behaviour — ^Log (^^^ at large u, 



(p+{uj,fi) = — — - 



/i^ 2{aB - 1) . f 
Log 



;iio) 



where the parameters A^^ and as are defined in terms of the integrals 



Ag^ = / du- 



UJ 



Jo UJ \ /i 



1111 



We realize that, due to the logarithm in the radiative tail, the moments for 
> are divergent. 



As explained below in Subsection 10.11, the parametrization (11 101) is a simplified 
and approximated form of a full QCD sum rules calculation that includes radiative 
corrections. 
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A thorough study of the radiative corrections has been done by Lee and Neubert 
[9]. The lowest positive moments M^^ and M^^ are given as a power expansion in 
the ultraviolet cut-off A[/y 



Mi°)(A, 



uv) 



1 + 



Cpas 



4:71 



16A Cpas 



3A 



uv 



An 



Log 



A 



uv 



-2 Log2 



'A 



uv 



/i 



+ 2 Log 



'A 



uv 



TT 

12 



112) 



uv) 



A 



4A f Ci^a^ 



-2Log 



-4 Log 



'A 



uv 



/i 



+ 6 



8 Log 



'A 



uv 



7 
4 



TT 

12 



;ii3) 



0, one recovers Mq = 1 and Mi = Lee and Neubert define the 



In the limit as 

radiative tail of the function (p^{uj,fi) by the prescription 

'dMo{Auv,l^y 



dA 



that gives 



UV 



4A 1 

6 UJ 



(114) 



2-Log ^ 



;ii5) 



Notice that for large enough u the radiative tail becomes negative, and that at 
lowest order it agrees with flllOp . 

To include the radiative tail is not without ambiguity. To add it to our calculation 
of the long distance part of (f^{uj), we follow two different models, that have different 
continuity properties, but that lead to almost identical results : first, a similar 
procedure to the one proposed by Lee et al. [9J, although different in its details; 
second, a procedure close to the one followed by Braun et al. [7]. We now expose 
both methods, and compare the results at the end. 



9.1 Model to add the radiative tail following Lee et al. 

Our first model for the function ip^i^uj), including the radiative corrections fol- 
lows essentially the prescription of [9] 



BT 



RAD, 



(116) 
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where ip^'^{uj) is the function (13T|) or (IMI) with the harmonic oscillator wave function 
( 1851) and the parameters of the Roman shape function ( |98l) . 

Without loss of generality at leading order in a^, we take the radiative tail 
ip^^^{uj) as given by expression (IllSp with A replaced now by Asf{i^), as determined 
from the fit to i? — Xg'j 



Cpas 



TCUJ 



+ 



4A5f(/u) 1 



Log ( - 



;ii7) 



118) 



;ii9) 



(120) 



The relation between Uf, defined by the vanishing of the radiative tail to ensure 
continuity, and AspifJ') is given by 

4A5f(/^) l_ " 
3 ut 



The relation between AsrifJ') = A5i?(/i,/i) and A is [25] 

A = A,H/^) - ^ 4/. 
We choose from now on /i = 1.5 GeV as an illustration, 
/i = 1.5 GeV Cpas = 0.470 

For this value of fi one has, from B —>■ X^j and B —>■ X^iVe [9l [33] , 

A5^(;x) = (0.65 ± 0.06) GeV (/i = 1.5 GeV) 
that is in agreement with the determination (1961) from 



log ( — 







1211 



■^5f(/^) becomes a function of Ut and is a parameter to be determined by the 
matching with the QCD behaviour ffTT2l) . ffTT3ll . 

In terms of A5^(/i), at the lowest order in one can rewrite the moments (I112p 
and ffml) 



1 + 



16A5i7'(/i) CjT^a^ 



3A 



uv 



An 



Log 



-2 Log^ 

Auv \ _ 
fx i 



'A 



uv 



fx 



+ 2 Loe 



A 



uv 



fx 



71 

12 



;i22) 



M^^\Auv) 



A 



uv 



4:71 



4:71 



-4 Log 
-2Log2 ( 



A 



uv 



Auv ' 



+ 6 
h 8 Log 



'A 



uv 



fx 



4 12 



CpcXs 16/i 

'^TT 3~ 



(123) 
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Notice that a new term proportional to as appears in the first moment m!^^ {-^uv) ■ 
Using the LCDA ^p^{uj) from the BT model fl34|) with the harmonic oscillator 
wave function (l85l) and the parameters (!98l) we then compute the moments in the 
model 



M. 



{O)model 



M 



{l)model 



(A 
(A 



uv) 



uv) 



did ^^{uj) 



(124) 



and match with the OPE expressions fll22p and fll23p . 

Making the approximation, that will be discussed below, 

duj if^iuj) ^ / du (p^{uj) 
Jo 

du uj '-pAoj) = / duj uj ^AoS) 
h 

the matching implies 



(125) 
(126) 



uv t 



r(0) 



(0), 



that gives, since M 



(1) 
+ 

(0)BT 



{\)BT 



(1), 



(127) 



M 



{\)BT 



N 



Equation (11280 gives N in terms of ut 

Cpa 



+ 



47r 

16AsF(/i) CpOis 



-2 Log^ p + 2 Log 



47r 



Log 



and equation (11291) gives, expanding to first order in 



M 



{1)BT _ 4A5i7'(/i) 



1 + 



4:71 



47r 



4 Log m - 6 



{-6 Log {^^ 

CpUs 16/i 



16A5F(/i) 



12 



Log 



(128) 
(129) 



(130) 



- 1 



(131) 
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Using the result fllOSp of the model M_\_' = 0.964 GeV, this equation gives 
another relation, besides (I12ip . relating ujf and AsrifJ')- 

From ( ]12ip and ( ]13ip we can solve for ut and AspifJ') for a given value of fi and 
see if the value obtained for AsrifJ') is consistent with the known value of Asrif^) 
(I120p from the fit to 5 — X^^y. Once Ut and AsrifJ') are known one can compute 
N and A from (fTSOj) and ffTTSD . 

From ffml) and ffml) we obtain, for /i = 1.5 GeV, 

ut = 3.288 GeV = 0.974 

Asf{i^) = 0.578 GeV A = 0.354 GeV (132) 

These values for Asrif^) and A that come from the OPE constraints are only about 
10 % lower from the values fll20p coming from the fit to i? — >■ Xg'j and B —>■ X^iVi. 
We conclude that the situation is good enough. 

A different status from A has the parameter A^^, that enters in a number of 
processes that we examine blow. The value that we obtain including the radiative 
corrections is 

A;^^ = 1.429 GeV"^ (133) 

to be compared with the value fll09p without the radiative tail. The correction is 
small. We find, for the parameter as ( lllip 

(Tb = 1.207 (134) 

The function flll6p with its radiative tail is plotted in Fig. 4. 



9.2 Model to add the radiative tail following Braun et al. 



This model follows the regularization of Braun et al. as illustrated by equation 
OllOp . but using the full radiative tail of Lee et al. (I115p . We set 



RAD I 



(135) 



where (p^'^{uj) is the same function (!3T]) or (IMP as in flll6p and we now regularize 
the radiative tail by making the replacement ^ T^Ti^i i-^- "we take 



71 Oj'^ + jJ? 
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. 6 



. 4 



0.2 



1 2 3 4 5 

Figure 4: The function ^p^{uj) in the BT approach with its radia- 
tive tail adopting the SF renormahzation scheme for fi = 1.5 GeV. 
The two curves show the two ways of gluing the radiative tail, follow- 
ing Lee et al. ^ (lower curve) or Braun et al. [7] (upper curve). 



+ 



4AsF(/i) 



UJ 



2- log 



-6 



(136) 



The radiative tails (11171) and (I136p differ at low uj. Now we do not have to determine 
the gluing point tUt. To impose the OPE conditions we proceed as follows. We have 
two conditions to fulfill 



Mr'^^'(A^y)=Mo°^^(A^v') 



M{"°'^^'(At;y) = Mf^^(A 



(137) 



These conditions give, respectively 
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(139) 
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log (- 
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where m|^^^^ = 0.964 GeV from ffTOSD . 

We have to solve these two equations for and A^^? in terms of Auy and find in 
which region of Auv these quantities become approximately constant. We find that 
this is the case for the reasonable range A[/y > 3 GeV, giving 



N = 0.856 



Asp{fi) = 0.563 GeV 



A = 0.338 GeV 



(140) 



i.e. a situation very close to the first model for AsrifJ')- 

Although the radiative tails are different, the final functions are almost identical, 
as shown in Fig. 4. The reason is that the positive part of the radiative tail in the 
latter model is compensated by the lower value of f ll40p . imposed by the OPE 
constraints. 

In this second model we have, for the quantities of interest. 



A;b^ = 1.432 GeV"^ 



as = 1.219 



(141) 



These values are very close to the ones of the first model (11331) . (I134p . 

We conclude that our results, due to the OPE constraints, are invariant relatively 
to the way of adding the radiative tail to the non-perturbative part. 



10 Proposals for LCDA in other approaches. 

The alternative theoretical method for the calculation of the LCDA is essentially 
the one of QCD Sum Rules. Here we will distinguish between the work that considers 
the LCDA at leading order and the one incorporating the radiative tail. Our aim is 
not a critical one, but only to show the great variety of ansatze that one can find in 
the literature for the functions ip^{u), and the corresponding varied results for the 
parameters A and A^^. 
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10.1 LCDA at leading order. 

A word of caution is in order here. In this part, that involves the LCDA at 
leading order, without radiative corrections, A is considered to be related to the 
center-of-gravity of the function Lp^{uj), i.e. to its first moment through 

4A , , 

<u>^ = — (142) 

Actually, in the calculation of the present paper, this first moment is assimilated to 
^^y^, since the parameters are obtained from the fit of the Roman shape function to 
B —>■ Xs7, that provides Asf- This will make easier the comparison with the other 
approaches in the absence of radiative corrections. 

Before going to specific theoretical schemes, it is worth to quote the bound found 
by Korchemsky, Pirjol and Yan [5J, that is independent of the precise form of ip^{uj), 
and follows from assuming positivity (p^{uj) > 

Ab' > i (143) 
^ - 4A ^ ^ 

Of course, the positivity condition ^p^{uj) > is violated by the radiative tail, as we 
have seen in the preceding section. 

10.1.1 QCD Sum Rules. 

Grozin and Neubert did obtain from the QCDSR result the simple form for the 
functions f±{uj) [2], 

^^^u;) = 4 e""/"" ^"^{uj) = \ e-"/"o (144) 

that satisfy relation (l39l) and give, for the positive moments (see also the analysis 
ofref. [31]), 

4A , , 

< u} >+ = 2 < uj = 2ujq = — (145) 

and for the parameter A^^ 

From the value A = 0.55 GeV used in this paper one gets 

Xb^ = 2.72 GeV-i (147) 
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On the other hand, Braun et ah [7] have obtained from QCDSR the following 
simple form for the normalized long distance shape of the LCDA ^p^{uj) 

V+{ij) = ^^u{2e,-uj)e{2e,-uj) (148) 

where Sc is the continuum threshold, and obtain 

4A , 1 3 9 , , 

<.>,= 2<.>_ = - A,' = - = ^ (149) 

that gives, for = (0.9 ± 0.1) GeV [7] 

A ^ (0.67±0.07) GeV A;b^ ^ (1.68 ± 0.18) GeV"^ (150) 

Also, within the QCDSR approach. Ball and Kou ^ find the above relation and 
the numerical value, adopting the value A = 0.68 GeV [351 ES] 

Ao^ = -1 ^ 1.67 GeV^^ (151) 
^ 8A ^ ' 

Both refs. [H [8] find the same expression (11481) . However, we must emphasize 
that in the work of Braun et al. [7], formula (I148p represents just the simplest 
contribution to the QCD Sum Rules. The full result of is much more complicated, 
and the formula (11 101) represents a simple parametrization of the full expression. 
In ref. [7] the long distance part and the radiative tail follow together from the 
calculation, and the latter is not added by hand. 

10.1.2 QCD Factorization models. 

Within the QCD factorization approach of Beneke et al. [3] for the calculation 
of charmless non-leptonic B decays {B Ktt, B tttt, ■ ■ ■) the following range is 
adopted : 

A;^^ = (3.5 ± 1.5) GeV"^ (152) 

This is a guess essentially based on the determination (I146P of [2J, using A = 0.4 GeV. 
The same range is adopted in [6j for the description of the decay B^ itu^. 



32 



10.1.3 pQCD Factorization models. 

A number of models have been proposed for the function $+(0 defined in ([3]) 
within the framework of the Perturbative QCD Factorization (pQCD Factorization) 
approach, needed in the description of non-leptonic charmless B decays. Integrating 
over p± or equivalently taking impact parameter b = 0, the following models were 
proposed [361 [STJ [38] ^oi the function $+(0 defined by ([3]), 

$^(0=iVF(Oexp 



1 I niB^ 



with a range of values for cj^ in the interval 



(153) 



0.25 <ujb< 0.65 



(154) 



and different models for the function F{^), 



F{i) = i\l-i)\ iil-i) or Ve(l-0 



(155) 



with N determined by the normalization condition (jl]). 

To compare with the present work we have to make the change of variables ([5]) 
and take the limit uih oo at fixed. Using relation ([7]) one finds the following 
models for ip^{(jj), 



i^l{uj) = N G{uj)exp 



2 \ujbJ 



(156) 



with 



G{uj) = LU , iu or \fu) (157) 

and determined by the normalization condition ([8]). 

Just a few comments are in order to compare with the present work. First, in 
the BT scheme with an harmonic oscillator potential, is not proportional to 

a gaussian, but given by the function fll06p . Second, the models fll53p - (11571) give 
generically low values for the first positive moment and very large values for the first 
inverse moment 



0.38 GeV < < >_ 



4A 



2.66 GeV"^ < < 



< 0.48 GeV 
A^^ < 4.18 



(158) 
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10.1.4 Results in the BT quark model. 

The Bakamjian-Thomas quark model of the present paper with harmonic oscil- 
lator wave function and parameters taken from the fit to the B — > Xgj spectrum 
using the Roman Shape Function yields 

<uj>+ = 2<uj>_ = = 0.964 GeV (159) 

and 

A;b^ = < uj-^ >+ = 1.521 GeV"^ (160) 

that satisfies the bound fll42p . The comparison between A^^ and < u is different 
than in the Lee et al. model (11441) . since we find approximately 

3 

< cj-^ >+ ^ - < (161) 

\~^^ is smaller than in the Lee-Neubert model due to the effect of the dynamical light 
quark mass 1712, that has a large value from the fit with the Roman Shape Function, 
and depresses < >+ due to the behaviour fH21) . 

From the value for Asf = 0.723 GeV from (11591) . not inconsistent with (I120p we 
obtain, for /x = 1.5 GeV, from (ITT9D and (fTTSD : 

A = 0.49 GeV (162) 

We summarize the values obtained for the parameters A and A^^, guessed or 
used in other approaches, that we compare with the results of the present paper 
(Fig. 1), in Table 1 and in Fig. 5. 

A remark is in order here in the comparison in Fig. 5 between ^p^{uj) from, e.g. 
the QCDSR result (I144p and our model. While in the former < u >+ = ^ with 
A = 0.55 GeV, we have in our model ([I59]), < to >+ = with Asf = 0.723 GeV, 
since it is Asf that is determined by the Shape Function. This explains why our 
(f^{uj) is more spread than (11440 . 
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0.5 1 1.5 2 2.5 3 



Figure 5: Different models for the function ip^{uj) in leading order. From higher to 
lower curves : the heavy quark limit of the three models of pQCD [36l [371 ISJ (]156p . 
ffTFTD with ub = 0.4 ; follows the models of QCDSR [2J with A = 0.55 GeV and 
[71 [8j (11481) ; the wider curve is the model of the present paper fl3T|) with the harmonic 
oscillator wave function fl85l) and the parameters fl98|) . that give A = 0.49 GeV. 



35 



Theoretical frame 


A 




Positivity 

ijD^(uj) > 15 




> 4 

— 4A 


QCDSR [2] 


0.55 GeV 


^ = 2.72 GeV"^ 

2A 


QCD Factorization 

PE] 


~ 0.45 GeV 


-2= = (3 5 ± 1 5) GeV"^ 


QCDSR [ZllH] 


(0.67 ±0.07) GeV 


4 = (1.7 ±0.2) GeV-i 

8A ^ ^ 


pQCD [36J-[38J 


(0.32 ± 0.04) GeV 


(3.42 ±0.76) GeV^^ 


BT model with 
Roman shape 
function parameters 


0.49 GeV 


1.52 GeV-i 



Table 1 : Results for the parameters A and A^^ in the different theoretical approaches 
for (p^{u}) in the absence of the radiative tail. 

10.2 LCDA with a radiative tail. 

Let us now discuss the impact of adding a radiative tail to the LCDA in the 
various approaches. Braun et al [7J have proposed the parametrization f lllOp of 
(p^{uj,fi) for /i = 1 including the radiative tail. Lee and Neubert [9] take as a model 
the long distance piece expression fll44p adding the tail (11171) with AsrifJ') replaced 
by AnAifJ-,!^), 

The consistency of (11631) with the first moments with the OPE (lll2l) -( fTT3ll im- 
poses constraints on and c^o that only depend on Ut and Ada = ^DA{fJ',fJ^) [9]. 
For /i = 1, using the relations between A/5^(yU, /i), AspifJ'*, f^*) and A one finds from 
A5F(/i*,/i*) = 0.65 GeV ([1201), Ada{i^,iA = 0-52 GeV (/i = 1 GeV), cut = 2.33 GeV, 

= 0.963, ujq = 0.438 and the results of Table 5. On the other hand, Lee et al. 
have shown that the model of Braun et al. for (f^{uj), for /i = 1 GeV is quite close 
to their own. 

We show in Fig. 6 and Table 2 the results of our models for adding the radiative 
tail of sections 9.1 and 9.2, compared to the results of [7] and [9]. 
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Figure 6: The function ^p^{uj) including the radiative tail for yU = 1.5 GeV. From 
higher to lower curves : model of the present paper ; Lee-Neubert model [9]. 



Method 


A (GeV) 


Ab' (GeV-i) 


Braun et al. [7j 


0.4- 0.5 


1.98 ±0.52 


Lee-Neubert |9J 


0.34 


1.86 ±0.17 


BT model (HD 


0.35 


1.43 


BT model (E]) 


0.34 


1.43 



Table 2 : Results for the parameters A and for /i = 1.5 GeV in the QCDSR 
approach and in our models ([1]) and ([2]) of sections 9.1 and 9.2 including the radiative 
tail. 

Some comments are in order here on the row of ref. [7J in Table 5. The range 
for A is the choice given in |7j, below formula (JT5l) . The value of A^^ (1.5 GeV) is 
obtained from X~^^ (1. GeV) (formula (39) of [7|) using {p^{u,fi = 1 GeV) (formulas 
(|39|) . fH3l) and fjlil) ) and the scale dependence for A^^(yu) given by fHT]) . 

11 Phenomena sensitive to the LCD A. 

Let us now review the observables that are related to the LCDA. In this respect 
there is essentially, the work based on QCD Factorization and the one based on 
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QCD Sum Rules. We will briefly review : 

(1) The decay B" -fiv^. 

(2) Hard scattering contribution to non-leptonic decay amplitudes like A{Bd —>■ 
TT+TT^) in the framework of QCD Factorization. 

(3) Asymptotic behaviour of the Isgur-Wise function C,{w) and the subleading 
form factor ^3(10). 

(4) Heavy-to-light form factors like i? ^ vr at = 0. 

11.1 The decay B" -fiVi. 

This decay is described by two form factors Fv{E^) and FaIE^) that, at tree 
level, are related, in the heavy quark limit, to the LCDA ^p^{uj) by 

Fv{E^) = Fa{E^) = <uj >+ = — (164) 

This process is directly related to the parameter A^^ and would be the most direct 
way of measuring it. Conversely, having a good theoretical estimate of A^^, the 
process '^tPi allows to measure \Vub\- 

A considerable effort has been devoted to the study of this decay going beyond 
the tree result (11641) . Korchemsky et al. [5] have computed the form factors for 
photon energies larger than Aqcd combining QCD methods for exclusive processes 
with HQET. They have written the leading twist form factors as the convolution of 
the B meson light-cone amplitude ^p^{uj) with a hard scattering term, and computed 
also Sudakov contributions. In a later paper, Braun et al. [7] have considered the 
radiative tail for the LCDA (f^{uj), reviewed in Sections 9 and 10. Descotes-Genon 
and Sachrajda [6] have studied the decay B^ 'yii'e in the framework of QCD 
factorization, demonstrating that indeed at the one loop order the amplitude can 
be written as a convolution of a perturbatively calculable hard-scattering amplitude 
with {p^{u). For the parameter A^, they use the guess of [3] A^ = (350 ± 150) MeV. 

The scheme of the present paper predicts a value A^^ = 1.43 GeV^^ that is 
in the lower range given in the different schemes of the literature, as reviewed in 
Section 10. This feature is due to the rather large value of the dynamical mass of 
the constituent light quark. 
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It is important to underline that the decay B~ — * 'yii'e is the only process 
that could allow, in principle, to directly measure A^^, modulo radiative and l/mg 
corrections. 

Recently, in a search for the decay 7^'*'7^, with i = e or fi BaBar has 

found the following upper bounds [39] , depending on the way of analyzing the data : 

Xb^ < 1.49 GeV"^ (165) 

or 

A;^^ < 1.69 GeV"^ (166) 

These bounds are fulfilled by our value (fT33D A^^ = 1.43 GeV"^ while they seem 
at odds with the predictions of other schemes (Tables 4 and 5). The predictions of 
refs. [7] and [9] are within Icr or 2cr in agreement with the bounds fll65p or f ll66p . 

11.2 Hard scattering in non-leptonic two-body B decays. 

The correction to factorization to the decays with an emitted light meson B — > 
D-K^-R-K, ... , that comes from a gluon attached to the spectator quark, called the 
Hard Scattering Amplitude, depends directly on $:^(0 or ^^{oj) [3|. In the case 
of the decays to two light mesons, it scales in terms of dimensionfull quantities like 
Gpc^s'm'b-^QCDy having the same behaviour as the leading term Gp^bJ^QCD- ^^is 
case of decays of B to two light mesons, this results in a contribution to the effective 
QCD factors like ai [3] 



C2 Cpnas , . 

ai.ll = ^ (167) 



with 



< (1 - x)-' >K < (1 - y)-' +rl< X'' >K XI 

(168) 

where is the well-known chiral enhancement factor [3]. This contribution to non- 
leptonic decays is proportional to A^^. BBNS [3J propose the number A^ = (350 ± 
150) MeV, based on the relation to A obtained in [2| and the bound [5J. The range 
adopted by BBNS A^^ = (3. 5 ±1.5) GeV~^ is larger than the prediction of the model 
of the present paper, A^^ = 1.43 or A^^ = 1.52 with or without the radiative tail. 
However, in the phenomenological analysis of two-body non-leptonic decays, as we 
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can see in expression fll68p . A^^ is afffected by a subleading although important term, 
chirally enhanced and proportional to the unknown logarithmically divergent factor 



lo T^^p{y)i where ^p{y) = 1 is the twist-3 pion light-cone amplitude. The second 
term in fll68p is parametrized phenomenologically by X'^ = (l + pj|/e*'^^)Log(mB//i) 
{Ph < 1), and fitted to the data. The modulus of the ratio between the two terms 
in ([TSHD is roughly |r-X^| = -^^i— ^|1 + p^e*'^«|Log(mB//i), of 0(1). This 
term results in a large uncertainty on Hktt (Fig. 5 of [3j), that for the real part is 
0.5 < KeHKn < 2.5, with a somewhat smaller uncertainty for the imaginary part. 
Therefore, strictly speaking, due to this unknown term, that plays a non-negligible 
role in the description of the data, two-body non-leptonic decays do not allow to 
make a model-independent extraction of the parameter A^^. We can conclude that 
data on non-leptonic decays into two light mesons, although dependent on A^^, are 
not a model independent determination of this quantity. 

11.3 Asymptotic behaviour of the functions ^{w) and ^'s{w). 

Grozin and Neubert [2j have given the behaviour of the elastic Isgur Wise 
function ^(cosh6') in the large recoil limit w ^ 1 (6' ^ 0), 

C 

e(cosh0) = levra,-^ f < lj-^ >+ < cj"^ >_ e"^^ (169) 



where 2/ = ju^'^Q- In our notation, this writes 

i{w) = vr (/bv^)' < a;-2 >+ < uj'^ >_ — (170) 

where 



, f°°duj p., . 1 duo p., . 

<uj-^>+= ^<^+(c^) <cj-i>_=/ — v^^(c^) (171) 

JO UJ'^ Jo 00 

From the description of the behaviour of the functions ^p^{uj) in the different limits 
in Sections 6 and 10, one realizes that these two integrals diverge for massless light 
quarks. However, the non-vanishing light quark mass of the BT model, implying 
the behaviour 0421) . provides a natural infrared cut-off. 

Using expressions (I3T1) or for non-vanishing light quark mass and the wave 
function fl85|l with the parameters fl98l) . one finds the finite results 

< >+ = 3.826 GeV-^ < uj-^ >„ = 2.322 GeV"^ (172) 
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From the value fs^/T^B = 0.388 GeV^^^ found in Section 8.2 in the present case, 
one finds the behaviour 

i{w) ^ 1.337 TT ^ (173) 



This gives, for /i = 1 {Cpcts = 0.624), an order of magnitude ^iw) ~ smaller 
but of the same order of magnitude as with a function i{w) = i^-:;;^^ and a slope 
of the order = 1. 

A comment is in order here. Grozin and Neubert [2] argue that the logarithmic 
singularities of < u^"^ and < > would be cut off by the transverse momenta 
and virtualities of the light quarks in the mesons. In our scheme, the transverse 
momenta do not play such a role, since the divergence remains taking those into 
account, as shown by the explicit expressions for ^p^{uj) in the massless limit that 
can be read from fl3T]) or fl34l) . It is the non- vanishing dynamical light quark mass, 
that one can consider as a "virtuality" of the massless quark due to the < qq > 
condensate, that ensures the finiteness of the moments < cu"^ >+ and < cu"^ >_. 

The subleading function ^^{w) [IQ], coming from l/mq perturbations to the 
current, behaves [2] in the large recoil limit 6^0 like 

U^oshd) = Anas ^ f < >l (174) 

or, in our notation, 

U^) = I ^ UbV^b? < u.-' >l ^ (175) 
and from the value (11331) or (11411) for < uj^^ >+ = A^^, 

< uj-^ ^ 1.521 GeV"^ (176) 

we find 

e3(^) = 0.174 TT^iGeV (177) 

that gives 

Uw) ~ ^ (178) 
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11.4 B ^ TT form factor at = 0. 

From eq. (l23l) of [1] one has up to terms in 1/(1 — u), 

f--(O) . / .5 ,u (179) 



''6 

i.e., 



' l\c Jo UJ Jo [I — UY 

However, this expression is not usable due to the infrared divergence of the integral 
over M, since ^t,{u) = 6m(1 — u). 

One way out of this problem is the use of the Light Cone QCD Sum Rules 
approach by Khodjamirian et al. that provides a simple and explicit expression 
for the heavy-to-light form factors in terms of LCDA. Taking as an example the 
form factor F^~^'^{Q) one has, 

^r'^(O) = 7^ r^^exp (--^) (^) (181) 

An approximation to this expression has been used in the first reference of |4], 

^'=,,exp(-^),?(-i-) ..«,0,f ..exp(-^) (182) 

that gives, from the integral relation (HOj) . valid only for vanishing light quark mass, 



exp 



^5 



M2 



;i83) 



Some remarks are in order here : 

(1) If (p^{uj) and ip^{uj) are related by (l40l) . the numerical results for F^~~*'"{0) 
from (11811) and (I183p are very close. 

(2) As pointed out in [4j and we observe from (I183p . for a given window of the 
Borel parameter and a value of Sq, 

0.5 GeV^ < < 1.2 GeV^ ^ = 0.7 GeV^ (184) 

F^^'^(O) is very sensitive to the precise value of X^^. The value adopted in [4] is 
A;b^ = 2 GeV-i. 

(3) The differential relation (HOl) between ^p^{uj) and ^p^{uj) is only valid for 
m2 = 0, and is badly violated in the BT model examined here, where the value of 
m2 extracted from the B X^'-f spectrum is rather large. 
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(4) A non- vanishing dynamical light quark mass m2 has a dramatic impact since 
then ip^{0) = 0. We are not allowed then to use (11831) but must rely on the relation 
(11811) . Using this expression and the parameters (11841) one obtains a very small 
value for F^^'"{0) ~ 0.015, the reason being that ^p^{uj) vanishes at a; = 0. One 
would need a much larger value of Sq to get an appreciable contribution of y:>^{uj) 
to the integral (11811) . Our conclusion is that with our prediction for f^^u) (Fig. 2), 
relation (11811) and the set of parameters (I184p . one cannot describe the form factor 
pB^-^ ^ This feature deserves further investigation. 

12 Conclusions. 

In conclusion, within the Bakamjian-Thomas relativistic quark model, that in 
the heavy quark limit yields covariant form factors and Isgur-Wise scaling, we have 
computed the B meson Light Cone Distribution Amplitudes ip^{(jj), that are also 
covariant in this scheme, and satisfy, in the limit of vanishing dynamical light quark 
mass, the integral relation given by QCD in the valence quark-antiquark sector. We 
have also computed the Shape Function S{uj) that enters in the description of the 
decay B — > Xs7. The Light Cone Distribution Amplitudes and the Shape Function 
are related in the BT class of models and given in terms of the Qq internal wave 
function. Using a gaussian wave function, we have shown that the shape function 
is identical to the so-called Roman Shape Function. Using the parameters of the 
latter that fit the B Xg'-y spectrum, we have predicted the LCDA (p^{u). We 
have discussed the role played by the dynamical mass of the light constituent quark 
and included the short distance behaviour of QCD for ip^^u). Compared to most 
schemes in the literature, our model predicts a rather small value for the parameter 
A]j^ = 1.43 GeV~^, due to the rather large value of the constituent light quark mass, 
fitted from the B Xs'j spectrum. This value for A^^ fulfills the upper bounds 
obtained by BaBar from the search of B^ 'yi'^i^e- Moreover, the non-vanishing 
constituent light quark mass has the important implication (p^{0) = 0. We have 
compared with other theoretical approaches and discussed the phenomena that are 
sensitive to the LCDA. 
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